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Abstract. In the present article we obtain the large A*' asymptotics of the par- 
tition function Zn of the six-vertex model with domain wall boundary conditions 
on the critical line between the disordered and antiferroelectric phases. Using the 
weights a = 1 — X, 6 = 1 + X, c = 2, |x] < 1, we prove that, as A ^> oo, Zn = 
CF^ N^/^'^ (l + 0(A~^)), where F is given by an explicit expression in x and the x- 
dependency in C is determined. This result reproduces and improves the one given in 
'~{ ' the physics literature by Bogoliubov, Kitaev and Zvonarev [7j . Furthermore, we prove 

that the free energy exhibits an infinite order phase transition between the disordered 
and antiferroelectric phases. Our proofs are based on the large A asymptotics for the 
underlying orthogonal polynomials which involve a non-analytical weight function, the 
Deift-Zhou nonlinear steepest descent method to the corresponding Riemann-Hilbert 
problem, and the Toda equation for the tau-function. 

1. Introduction and statement of the main result 

We begin with the description of the model under consideration. Given a square 
lattice in 1? with N"^ vertices, we assign arrows along each edge obeying the following 
rule: At every vertex two arrows point in and two arrows point out. Such a rule is called 
the ice-rule and it only admits six possible arrow configurations, see Figured! hence the 
name of the model: the six-vertex model or the model of two- dimensional ice. 

On the boundary of the lattice we impose the domain wall boundary conditions 
(DWBC), that is all arrows on the top and bottom side of the lattice are directed 
inside the lattice and all arrows on the left and right side point outside, see Figure [2] for 
a possible arrow configuration with DWBC. 

To each of the six possible arrow configurations {i) we assign a positive real-valued 
Boltzmann weight tt;^, i = 1, . . . , 6 and define the partition function of the model as 

configurations i=l 

with Ui denoting the number of vertices with arrow configuration (i). By definition, 
depends on six parameters: the weights Wi. However recalling some conservation laws 
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Figure 1. The six allowed arrow configurations at a vertex 



Figure 2. An example of a possible 4x4 configuration with DWBC 
(see for instance [I], [1] or [12]) we can reduce the number of parameters to two, namely, 



ZNia,a,b,b, c, c). 



(1.1) 



where 

and furthermore 



a = ^WlW2, b = y/W3W4, C = y/W5WQ , 

Zj\[{a, a, b, b, c, c) = Zn 



I a a b b 

c c c c 



(1.2) 
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Thus, a general weight reduces to the two parameters -, -. In order to study the phase 
diagram (see Figure [3]) of the model, we introduce the parameter 

2ab 

We distinguish the following three physical phase regions: the antiferroelectric phase 




1 



Figure 3. The phase diagrom of the six- vertex model, with (F), (AF) 
and (D) denoting the relevant phases 

(AF), A < —1; the disordered phase (D), — 1 < A < 1; and, the ferroelectric phase (F), 
A > 1. In all phases the weights are usually parameterized in the following way: For 
the ferroelectric phase, 

(F) o = sinh(t-7), 6 = sinh(t + 7), c = sinh(27), < |7| < t, (1.3) 
for the antiferroelectric phase, 

(AF) a = sinh(7-t), 6 = sinh(7 + i), c = sinh(27), < |t| < 7 (1.4) 
and for the disordered phase, 

(D) a = sin(7-t), 6 = sin(7 + t), c = sin(27), < |t| < 7 < ^ . (1.5) 

Here we discuss the critical line between the disordered and antiferroelectric phase, hence 

A = -1 

which corresponds to 

a b 

- + - = 1. 

c c 

Instead of using the disordered phase parameterization p.Sp in the limit 7 — t- 0"*", we 
choose the following parameterization for the weights a, b and c: 

a = l-x, b = l + x, c = 2, xG(-1,1). (1.6) 



EXACT SOLUTION OF THE SIX- VERTEX MODEL 



4 



The above described six-vertex model with DWBC was introduced by Korepin in [15] 
and further studied in [H] and [8]. Those works lead to a remarkable determinantal 
formula for the partition function with DWBC, which, with weights parameterized as in 
dm), is 



TN, (1.7) 



where tn is the Hankel determinant, 

/ \N c 2 

The determinantal formula implies that r^v solves the Toda equation, 

2 d 
tnt'n - (r^) = TN+iTN-i, n > 1, (') = —. (1.9) 



As was noticed by Zinn- Justin [19], the Hankel determinant r^r can be expressed in 
terms of the partition function of a random matrix model, and then it can be reduced 
to orthogonal polynomials. On the critical line under consideration, the expression in 
terms of orthogonal polynomials can be obtained as follows. We write f{x) as the 
Laplace transform of a continuous measure, 

oo 

^^"^^ = {l-x){l + x) = I ^'■''"(^)'^^' "^(^) = 

— oo 

and define the potential V{X) = |A| — xX. This implies for tat 

e- EL. ^(A.) f Yl A|-M det {^~%^^X, • • • dA^ 

— oo — oo ' ^ 

and thus for any permutation a G Sn acting on {Xi}^^ 

e- EL. nA.) f Yl sgn{a)Xl-^l^ ) det (Ar^);;;.^,^Ai • • • dA^, 

— oo — oo ' ^ 

hence after summation and identification of a Vandermonde determinant 

oo oo 

^^^m I "' j ^~^'='^^'''^T{i^^-^3fdXi---dXN. (1.10) 

— oo — oo ^^-J 

Up to the factor ^^^7, the expression on the right is the partition function of eigenvalues 
in the random matrix model with the potential V . 

Introduce now monic orthogonal polynomials {pn(s)}n>o with respect to the measure 
d//(s) = e~^'^^^ds on the real line, 

00 

j Pnis)Pm{s)dfl{s) = K5nm ■ 
—00 





EXACT SOLUTION OF THE SIX- VERTEX MODEL 



5 



Then, by using the orthogonaUty condition, one can simphfy (jl.lOp to 

N-l 

rN=ll hk. (1.11) 

fe=0 

Thus, for the partition function, via (|1.7|) . we obtain the formula 

N-l , 

= («^)"^ n (1-12) 

We will prove the following asymptotics of the normalizing constants /i^v: 
Theorem 1. As N ^ oo, 

^ " ^ 'i + t4+-^)> (1-13) 



(iV!)2 \^2cosf^7 V 12N 



(-l)^cos(^x(iV+ i)) r 21nlniV 1 - 27 - 41n2 - 21n (cos f^) 
2iV(lniV)2 IniV hT/V 

|+0(iV-2), (1.14) 



InlniV^ ^ 



where 
En 



InN 

and the error terms in (jl.l4p are uniform on any compact subset of the set 

{x G M : -1 < X < 1}. (1.15) 

The proof of Theorem[T]is based on the Riemann-Hilbert approach for the potential V 
with a singularity at the origin. V is close to the class of the Freud potentials considered 
in the paper of Kricherbauer and McLaughlin [T7] . Our approach is somewhat different 
from [l7], and it gives a more detailed asymptotic formula for hf^f, which is needed in 
the asymptotics of Z]\f. Our second main result in the present paper is the asymptotics 
of Zj\f. 

Theorem 2. On the critical line between the disordered and antiferroelectric phase re- 
gion with — l<x<l, as N ^ 00, 

Zn = CF^'nt2 [1 + 0{N-^)) , (1.16) 

where 

1 

r.2 



7r(l - X^) f 



^=^^- C^Coj^cos-j , (1.17) 

and Cq > is a universal constant. The error term in (|1.16p is uniform on any compact 
subset of the set (jl.lSp . 

Remark. The leading term in formula (|1.16p for Z]\f, after a rescaling, coincides with 
the one obtained in the physics literature by Bogoliubov, Kitaev and Zvonarev [?]■ The 
error term estimate of 0{N~^) is much stronger than a logarithmic estimate which can 
be derived from [7j- Expression p.l7p for the constant C is completely new. 

The current article is a continuation of the work of the first author with Vladimir 
Fokin in [3] and with Karl Liechty in [3], [5] and [6]. They prove conjectures of Paul 
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Zinn-Justin in [19] on the large N asymptotics of in different phase regions: In the 
disordered phase with parameterization (jl.Sp for some e > 0, 

(D) Zn = CF^^N''(l + 0{N-^)y N^oo (1.18) 

with 

7r(sin(7-t)sin(7 + ^)) 

^ 12 37r(7r-27)' 27 cos ^ ' ^ 

and an N independent constant C > 0, whose value is unknown. In the antiferroelectric 
phase with (jl.4p 

(AF) Zn = CF^^^4{Nu)(l + 0{N~^)y N ^ oo, (1.20) 

where 

vr/ t\ 7r(sinh(7-t)sinh(7 + t))i?;(0) 

2 V 7/ 27t?i(tj) 
and t?i(2;) = i?i(z|(7), i?4(z) = '!94(-z|q') are the Jacobi theta functions with the elliptic 

nome q = e ^7 . The constant C > is also here unknown. 
Furthermore, in the ferroelectric phase ()1.3p for any e > 

(F) Zjv = CF^'G^(l + 0(e^^"^)), N ^ oo 

with 

C = l- e-^^, G = e^-*, F = sinh(7 + t). 

Finally on the critical line between the ferroelectric and disordered phases with the 
parameterizations a = ^(t — 1), 6 = ^(t + 1) and c = 1 with t > 1 

(F-D) Zjv = CF^'G^iV3(l + 0(Af-i/2)), iV 

where 



oo 



G = exp 



with C(s) denoting the Riemann zeta function and a positive, yet unknown, constant 
C > 0. 

To calculate the order of the phase transition on the critical line between the disordered 
and antiferroelectric phase regions, we would like to compare the free energy in these 
phase regions as we approach the critical line. First we have to adjust the partition 
function in the regions. Observe that as we approach the critical line, the parameters 
t,7 — ;> 0, hence c ^ 0. We rescale Zj^ according to (|1.2p : 



f a a b b \ _j^2 
Zn 1, 1 = c Z^ia, a, b, b, c, c). 

\c c c c J 

In the disordered phase we obtain from (|1.18p that 

Z;vf-, 1, l) =CF^'iV''(l + 0(iV-)), iV^oo, (1.22) 
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with the same C and k as in (jl.lSp . and with 



Here a, b, c are parameterized as in (jl.Sp . 

In the antiferroelectric region we obtain from ()1.20p that 

Zm(-, =CF^'MNoo)(^ + 0{N~')), N^oo, (1.24) 
\ c c c y ^ ^ 

with the same C and a; as in (|1.20p . and with 

(AF) F = F(y.t)=-^^^. 0<W<y. (1,25) 

Here a, b, c are parameterized as in p.4p . Let us remind that 

oo 

=i9i(z|g) = 2^(-l)"g('^+5)%in((2n + l)2), g = e" 



2l£ 

27 , 



n.=0 



1 + - 1. 



(1.26) 



7. 

Observe that as (7 — 0, 

??i(z|g) = 2gi/^sinz + 0(g9/^), ^[{0\q) = 2q^/^ + 0{q^/^), (1.27) 

hence 

tr^ = + = ^ + (1.28) 

■Wi^tj) smw cos 1^ ^ ^ 

The free energy in the antiferroelectric region can be written as 

F(7,t) = F-s(7,t)+F^'"S(7,t), (1.29) 

where 

27c cos (1.30) 
a = sinh(7 — t), 6 = sinh(7 + t), c = sinh(27), 

and 



27cV^?i(w) cosf^. 

Observe at this point, that F{j, t) in the disordered phase region see (jl.23p and -F'''^^(7, t) 
in the antiferroelectric one see p.30p are given by the same expression, but the underlying 
parameterizations of a, 6 and c are different. 

Introduce now the coordinates (x, y) on the phase diagram as 

a I — X b I + X /-.-.N /-,o-.^ 

- = ^^ + y, - = ^^ + y, xG(-i,i), (1.31) 

c 2 c 2 

so that y > corresponds to the disordered phase region, whereas y < to the antifer- 
roelectric one. Employing the change of variables 

(7,t)^ (x(7,i),y(7,i)), (1-32) 
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defined by (jl.3ip and (|1.4p in the antiferroelectric phase region, and, respectively, by 
(jl.Sip and (jl.Sp in the disordered phase region, we set 

FD{x,y) = F{-f{x,y),t{x,y)) in (D), FAF{x,y) = F{-f{x,y),t{x,y)) in (AF), 

and on the critical line (cf. Theorem [2]) , 

Also, we set in the (AF) region 

F^<|(x,2/) ^ F-g(7(x,y),t(x,y)), F^°S(x,y) ^ F^-s(7(x, y), t(x, y)) . (1.33) 

We prove the following result on the order of the phase transition on the critical line, 
see also the work 1191 of Zinn-Justin. 



Theorem 3. We have the following properties: 

(1) For any x € (—1, 1), change of coordinates p.32p implies that as y ^ 0, 



^ = 2 V (1 + 0{y)) , t = 4x^ ^ {l + 0{y)). (1.34) 

(2) The functions F£){x,y) and F^p(x,y) can be analytically continued in y to a 
neighborhood of the point y = 0, and in this neighborhood 

FD{x,y) = F\'^{x,y). (1.35) 

In addition, we have that 

Fd{0,x) = Fc{x) = F'^^{0,x). (1.36) 

(3) As y ^ 0, F2p^{x,y) is exponentially small, 



AF 

AF 



F^7(x,2/)=0(e--), XG(-1,1), (1.37) 



that is we have an infinite order phase transition. 

The setup for the remainder of the paper is as follows: We will first proof Theorem [T] 
within the framework of the Riemann-Hilbert approach to orthogonal polynomials. We 
start with a rescaling of the weight and the construction of the equilibrium measure. 
Afterwards a sequence of transformations is carried out which allow to approximate the 
global solution of the Riemann-Hilbert problem by local model functions, and enable 
us to solve the initial problem via iteration. After that, Theorem [2] is a corollary to 
Theorem [H via (jl.l2p and (jl.9p , with some extra arguments enabling us to prove the 
0{N~^) estimate of the error term. The proof of Theorem [3] will be given at the end of 
the article. 
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2. ReSCALING and THE EQUILIBRIUM MEASURE 

Let us scale the variables in (jl.lOp as Aj = A^/ii, so that 



oo oo 

/ ••• / e-^Sf.inM.)JJ(^^_^^.)2rf^^...rf^^. (2.1) 



Now use the rescaled monic orthogonal polynomials {pN,n{s)}n>o with 



thus satisfying 



and simplify (|2.ip similarly as in (jl.lip 

Af-l 

= n 

i=0 

Our original task of studying the asymptotics of orthogonal polynomials with respect to 
the weight function w{s) = e~^^^^ has therefore simply been transfered to polynomials 
with the underlying weight 

w{s) = e-^^'-'\s£R; V{s) = \s\ - xs, \x\ < 1. (2.2) 

Our strategy will now focus on the determination of the large N asymptotics of the 
normalizing constants /jatvi and then, via = N'^^'^^h^ N, the evaluation of tn from 

dm]). 

To this end introduce the normalized counting measure v on R 

1 ^ f 

and rewrite the integrand in ()2.ip as 

i<j 

with 

H{v) = j j \n\t- s\^^du{t)dv{s) + j V{s)du{s). 

Since we are particularly interested in the behavior of r^v as — )• cxd, we expect the 
value of Tat in (|2.ip to be focused in a neighborhood of the global minimum of the 
energy functional H{v), with v reaching over 

M'^iM) = {pe Borel measures on M : 
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In \t - s\^Ufi{t)di2{s) + / V{s)dn{s) 



= inf 

has a unique solution = ^ € 7W^(]R), called the equilibrium measure. 

In the given situation the underlying potential V{s) = \s\ — xs is convex. In this case 
the support of the equilibrium measure consists of one interval, 

J = supp(/i^) = [a, P] C M. (2.3) 

For our purposes the equilibrium measure will be essential in the construction of the 
g -function 



g{z) = / ln{z - w)dir {w) = / ln(z - w)p{w)dw, z G C\(-oo, f3q], (2.4) 



where we choose the principal branch for the logarithm. Moreover the (^-function in turn 
determines the equilibrium measure uniquely by the following variational conditions: 
there exists a real constant I, the Lagrange multiplier, such that 

z£R\J: g+{z) + g^{z)-V{z)-l<0, z € J : g+{z) + g.{z) -V{z) - 1 = 0. (2.5) 

This characterization of the equilibrium measure leads to a Riemann-Hilbert problem 
for g'{z) which can be solved as 



.'(^) = ^/^^^, ^eCU (2.6) 



J 



It turns out that for the given potential, V{s) = \s\ — xs, the ^-function, its support, 
the density function p{z), and the Lagrange multiplier I can all be evaluated explicitly, 
see [T9] and p]. We have from ([23]), jZH) and ([22]) that 



,, , y^{z — a){z — /3) f sgn(w) — X dw 
9{z) = — 



2vri J y/(w - a){w - (3), w - z 

To evaluate this integral we use the residue theorem, 

dw 1 f diu 



yj{w — a){w — /3),{w — z) 27 y^(tt; — a)(zi; — I3){w — z) 

1 1 

iri vesw=z — , h ttz res< 



^{w — a){w — (3){w — z) -sj [w — a){w — I3){w — z) 

Tri 



a/(^-«)(^-/3) 
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where Cq denotes a closed contour around the interval [a, (3] such that both points z and 
oo lie to the left of Cq. Similarly 

^ 

sgn(?i;) dw ni ^ ^ . f dw 



yj{w - a){w - 13) w - z ^[z - a){z - /3) J {w - a){P - w){w - z) 
and the latter integrand has an exphcit antiderivative. We conclude 



= + 1 m - - ^\ .eC\|»,fl (2.7) 

2 vri ^z{(3 - a) 



choosing the principal branch for the square roots in (j2.7p . In order to evaluate the 
endpoints a, (3 of the support, use the resolvent 

P 

I, . f d/i^ (w) 1 ^, _2N 

g'{z) = / ^ =- + 0{z Z^OO 

z — w z ' 



and compare with the expansion in (|2.7|) as z — )■ oo. One obtains 

TT 

= vr tan —(1 + x) 

2^ (2-8) 



a = — TT tan — (1 — x), p = vr tan — (1 + 



(—a)/? = vr , /3 — a 

cos — 

For the ^-function itself we use an indefinite integration by parts, 
(7(z) = zg'{z) — J zg"{z)dz + C, C = const, 

where from ()2.7p 

/// N / Iz- P jz-a 



zvr(/3 — a) \\ z — a \ z — (3 
Thus, 



g{z) = zg'{z) + 2 In — a + Vz-/3) - (l + 21n2), zGC\[a,/3], (2.9) 
and the constant C = — (1 + 2 In 2) is obtained from a comparison of the asymptotics in 

g{z) = j ln{z- w)dfi^{w)= In z + 0{z~'^), z ^ oo (2.10) 

a 

with the asymptotics in ()2.9p . We move on to the density and use the resolvent equation 
again: 

g'4z) - g'_{z) = -27rip{z), ze[a,f3]. (2.11) 

However, from 



jl(~)-^g"("^-" , 2 V/3(. - g) + V-a(/3 - z) 
^ ^ 2 vri 
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Since 

V/3(z - a) - V-a(/3 - z) _ Vkl(^-a) 



Hence, via ()2.1ip . 

2 yPjz - g) + V-a (/3 - z) 



p(^) = £,ln^^^^^^-^^:^^:4^^> -e(a,/3), (2.12) 



thus the density has a logarithmic singularity at the origin. Finally, for the Lagrange 
multiplier, we state first, as a consequence of (|2.11|) and (j2.5|) . that for z G [a, /3], 

g+iz)-g.{z) = 2TTiJp{w)dw, g+iz) + g.iz) - Viz) - I = 0. (2.13) 
z 

But for such z 



-2ln{V z — a — i 13 — z) + C 



and 



9-iz) = — rln , + 21n [Vz - a - (3 - z) + C, 

hence 

/ = 21n(/3-g) -2(l + 21n2). (2.14) 
We have now gathered enough information to start the asymptotical analysis. 

3. RiEMANN-HlLBERT PROBLEM FOR ORTHOGONAL POLYNOMIALS 

The following Riemann-Hilbert problem (RHP), originally introduced by Fokas, Its 
and Kitaev in [13], is essential: Find a 2 x 2 piecewise analytic matrix- valued function 
Y{z) = Y'^"'\z) such that 

• y(")(z) is analytic for z G C\IR 

• Orienting the real line from left to right the following jump relation holds 

Yr\z)=Yl-\z)(l z^R 



• As z — )• oo, we have 

yW(z) = (/ + 0(z-i))z"'^3^ (73= 
The stated problem has a unique solution Y^^\z) given by 



1 

-1 



^ ' Wivrn-i(^) ^47r„_i(.)u;(.) '^^ 
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where 7r„(z) = + . . . is the n*^ monic orthogonal polynomial with respect to the 
measure dfj,{s) = 'w{s)ds and 

oo 



h 



n 



Furthermore Y^"'\z)z "'^^ admits a full asymptotic expansion near infinity of the fol- 
lowing form 

^(n) ^(n) 

and since in our situation w{s) = e~^^^^\ we obtain 

/i^,„ = -27ri(y/"))^2 (3-1) 

which, in terms of the previous discussion in the last section, shows that we need to solve 
the given RHP for Y{z) = Y^^\z). This solution will be obtained in the framework 
of the Deift-Zhou nonlinear steepest descent method [11] using techniques developed in 
|10| . These techniques allow to approximate the global solution Y{z) by local model 
functions, parametrices, and the iterative solution of a singular integral equation. We 
will elaborate the required steps in the following subsections. 

4. First transformation of the RHP - Normalization 
Recall (j2.10p and make the following normalizing substitution in the orginal Y-RHP 

Y{z)=exp(^^a3^T{z)exp(^N(^g{z)-^y3y z e C\R. (4.1) 

This leads to a RHP for the function T{z) 
• T(z) is analytic for z G 



The above properties of g{z), see (|2.13p . imply the following jumps 

/e-^(9+-9-) 1 \ 

T^iz) = T.{z)( ^ ,^(.,-.-)Me[a,/3] 



and 



/i N{g++g_-V-l)\ 

r+(z)=r_(z)( J " ^ ), zGR\[a,P]. 



• At infinity, T[z) is normalized 

Tiz) = I + 0{z^^), z^oo. 
Let us take a closer look at the listed jump matrices: For z S +00) 



g+{z) + g.{z) - V{z) - I = —aig^ ^ / ^ + 41n ^ ^ 
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and for z G (—00, a) 



g+{z) + g^{z) -V{z) -I = 2z H arg- ^ ^ / y j 



+4 In 



vr v^-z(/3-a) 



thus (see (IX^ ) 

5+(z) + (7-(^)-V^(^)-^<0, zGM\[a,/3] 
and we conclude for z G (—00, a — r/) U (/3 + 77, +00), r/ > fixed, 

'1 e^(^^++^^--^-') 



1 , ,/, iV^oo (4.2) 

where the stated convergence is, in fact, exponentially fast. Secondly the line segment 
[a, /3]: The given potential V{s) = \s\ — xs,s G M admits analytical continuation in the 
left and right halfplane separately via 

^^^)-\ -z{l + x), Rez<0, 

and this continuation is however two- valued on the imaginary axis. Since g~{z) = 
V{z) — g+{z) + l,z & [a, /3], the function 

G{z) = g+{z) - g.{z) = 2g+{z) - V{z) - I 

admits analytical continuation in a neigborhood of the segment [a, f3] into the first and 
second quadrant of the complex plane using the appropriate continuation of V{z) in 
(jOj) . On the other hand from (f233|) 

G{z) = 2-ni j p{w)dw (4.4) 
2 

and therefore 



ay y=Q vr ^J\z\{l3 - a) 

i.e. the stated local continuation of G{z) into the first and second quadrant satisfies 

Re G{z) > for Rez / 0, Imz > 0. (4.5) 

In the lower halfplane we argue in a similar fashion, indeed 

G{z) = -2g_{z) + V{z) + l 

admits local continuation into the third and fourth quadrant satisfying 

Re G{z) < for Rez / 0, Imz < 0. (4.6) 

These continuations will now be used in the following matrix factorizations 

^^~N(9Uz)-9-i-)) 1 \ _ f 1 ^\ f A f ^ 0^ 

= SljSpSl2, 
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motivating the second transformation of the RHP. 

5. Second Transformation of the RHP - Opening of lenses 

Let denote the upper (lower) lense, shown in Figure [H which is bounded by the 
contour j^. Define 



S{z) = <( T{z)Sl„ zgC^UC^, 
T{z), otherwise, 
then S{z) solves the following RHP 



(5.1) 







' £+ 


^ 7r 


a \ 

73" 




^4 











Figure 4. The second transformation - opening of lenses 

• S{z) is analytic for z G C\(M U F), with F = U 7^ U 73" U 7^: U (-ie, ie) and 
< e < 1 remains fixed. 

• The following jumps hold, with orientation fixed as in Figured] 

z G7+U72+, 



\^-NG{z) J, 

zG [a,0)U(0,/3], 



5+(z) = 5„(z)< 
Furthermore on the vertical line segment oriented upwards 



eiVG(-) ?)> ^ e 73 U74 



where 

with 

h{z) 



31 



S^{z) = S.{z)[ ^^^^^ Jj, zG(0,ie) 



{z) = e-^^-(^) - e-^^+(^) = -2fe-^-^i(^) sin(iViz) 



4z y^/3(z - a) + y^-Q(/? - z) ^ y/z- a - /3 - z 



2z^ 
+z mz 

TTl 



\//5 - a 



and similarly 



(5.2) 



(5.3) 
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with 



where 



4(-z) y^/3(z - q) + v^-a(/3 - z) ^z-a- - z 

J2{z) = —In + 41n 

m ^Jp — a — a 

_,_?(z£) ,„(_,). 

• At infinity, S{z) = I + 0[z~^) , z — )• oo. 

With ()4.5p and ()4.6p . the jump matrices will approach the identity matrix exponentially 
fast on the upper lense boundary 7^ and on the lower lense boundary 7-^ U7J . Also, 
as we have seen earlier in (|4.2|) . a similar statement holds along M\[q,/3]. Let us direct 
our attention therefore towards the vertical line segment (— ie, ie). For z = iy,0 < y < e 
we have 

ji{iy) = 2ie-^^-^y^''y+y^^^y^+^^^y^'^ sm{Ny) (5.4) 

where we introduced 



, . 4 y/P{iy - a) + y/-a{/3 - iy) yjiy-a + - iy 

hi{y) = -ln^ 7^=^ , /i2(y) = 41n 7^= 

TT -y/p — a VP — CK 

and with z = —iy, < y < e also 

j2{-iy) = 2ie-^''-^y^''y+y^'^y^+^^'-y^^ sm{Ny) (5.5) 

with 



^^^^ vr 
Hence, combining ()5.4p and ()5.5p 

ji(z) = ai(iy)e-^J^I'''^l sin(iVy), j2(z) = a2(-iy)e-^^l ^"(^^^l sin(iVy), < y < e 

with functions ai{z) which are analytic in a full neighborhood of the origin, satisfyting 
ai{0) 7^ 0. As a result of an optimization consideration, the function 

f{y) = e"-^y'^ysmiNy), yG[0,e) 
satisfies the important estimation 



/(y) = 0(^^J, N^oo, ye[0,e), 

hence also on the vertical line segment, the given jump matrices will, as A'^ — )• 00, 
eventually approach the identity matrix. All together we expect, and this will be justified 
rigorously, that as — t- 00, S{z) converges to a solution of the model RHP, in which 
we only have to deal with the constant jump matrix on the punctured line segment 
[a,0) U (0,/3]. Let us now consider this model RHP. 
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6. The Model RHP 

Find the piecewise analytic 2x2 matrix valued function M{z) such that 

• M{z) is analytic for z £ C\[a,/3] 

• Along [a,/3], the following jump holds 

M+(z) = M_(z)(„Vo)> ze[a,P] 

• M{z) has at most logarithmic singularities at the endpoints z = a, /3 

• M{z) = I + 0{z-^),z ^ oo 

A solution to this problem can be obtained by diagonalization 
with 

1/4 

5{z 



z — a 



defined on C\[q, (3] with its branch fixed by the condition (fE^)^^^ — t- 1 as z — >• oo. 

7. Construction of a parametrix at the edge point z = (3 
Fix a small neighborhood U of the point /5 and observe that 



^ Vz(/3 - a) 

\/z- a + a/z - /3 



+41n- 



as z S iY, z > /3. Simultaneously, 

G(z) = 2g{z) - V{z) - I = -co{z - pf'^ + 0[{z - (if), 

as z G n 7j*" and 

G{z) = -2g{z) + V{z) + / = co(z - pf'^ + 0{{z- fSf) , zGUHj^ 

where the function (z — /3)^/^ is defined for z G C\(— oo,/3] and fixed by the condition 

{z - Pf/^ > if z > /3. 

The stated local behaviors suggest to use the Airy function Ai((^) in the construction of 
an edge point parametrix. To this end first recall (see for instance [2]) that the function 
Ai(^) is a solution to the Airy equation 

w" = zw 

uniquely fixed by its asymptotics as C — ^ oo and — vr < arg ^ < vr 

Ai(C) = ^^e-l^"' f 1 - + ^^-6/2 + 0(^-9/2 

2^ \ 48^ 4608^ 
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as well as for — ^ < arg C < — ^ 



Ai(C) 



-1/4 



Now introduce for (" G C 



_2a3/2 
e 3'- 



48^ 4608^ V 



■63 



48^ 4608^ 



C — ^ oo. 



^o(C) 



AAi(C) e^^AAi e-^C 



Ai(C) e*3Ai(^e^*-C 
and observe from the previously stated asymptotics that as C — t- oo for < arg (" < vr 



(7.1) 



^o(C) 



0-3/4 



1 i 



1 



1 6i 

-1 



+ 



35 



4608(6/2 



-1 I2i 
-m -1 



+o(c 



-9/2^ 



On the other hand if — vr < arg ( < 0, we have 



MO 



2^ V 1 ^ 



1 



1 6i 



48(3/2 -V 4608(6/2 



+ 



35 



-1 I2i 

-m -1 



+o(c-^/2) 



e 3 



1 



We assemble the following model function 

r MO. 



\RH 



iO 



AoiO 

MO 



1 
-1 1 

1 -1 
1 



C — )■ 00. 



argCG (0,f ), 
arg C G (f ,vr), 



argCG (-f ,0), 



(7.2) 



MOC, -/ 



leading to the RHP depicted in Figure [5] 

• A^-^(C) is analytic for ( G C\{arg ( 



arg C G (-vr, 



.2zr Q 27r I 
3 , u, 3 , '1 ; 



27r> 
' 3 ' 



The following jumps hold, with the contours oriented as shown in Figure [S] 

27r 

arg C = =F^ 



arg C = 
arg ( = TT 



A^iO 


= A^^^iC) 


r 1^ 





A^^iC) 


= A^^iO 




:j 


AnO 


= A^^iC) 
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Figure 5. The model RHP near z = (3 which can be solved explicitly 
using Airy functions 



In order to determine the behavior of A^^ {(^) at infinity we make the following 
observations. First let arg ( E (— vr,— ^) and consider 



63 



1 

4^3/2 



1 



however here Re(C^/^) < 0, hence the given product approaches the identity 
exponentially fast as C — ^ co- Secondly for arg C G (t"' ^) ■ 



e 3 



1 

-1 1 



1 

4^3/2 
-63'' 1 



and also in this situation Re(^^/^) < 0. Both cases together with the previously 
stated asymptotics for ^o(C) imply therefore 



c 



0-3/4 



1 i 



1 + 



1 Qi 



+ 



35 



48C3/2 \Qi -I J 4608(6/2 \-l2i -1 



-1 I2i 



e-iC^/V3. 



(7.3) 



as C — 7- 00 in a full neighborhood of infinity. 

The model function {(,) will be useful in the construction of the parametrix to the 
solution of the ^-RHP in a neighborhood oi z = fi. We proceed in two steps. First 
define 



3iV\ 



2/3 



2g{z)+V{z) + l 



2/3 



\z- f3\ < r. 



(7.4) 



This change of variables is locally conformal, since 

2N ^ 



C(^) 



(- 



(z-/3)(l + 0(z -/?)), |z-/3|<r, 



and it enables us to define the right parametrix U{z) near z = /3 by the formula: 

U{z) = 5,(z)(-iV^)A^^(C(z))e§<^'^'(")"^ \z-P\<r (7.5) 
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with (^(z) as in ()7.4p and the matrix multipher 

-0-3/4 



-I I 



z — a 



Br{z)=[{ ; C(^)^ , Sr(/3)= / ;)(— (/3-a)) • (7.6) 



/3 



-<T3/6 



By construction, in particular since Br{z) is analytic in a neighborhood of z = /3, the 
parametrix U{z) has jumps along the curves depicted in Figure [6l and we can always 
locally match the latter curves with the jump curves of the original RHP. 

1 0^ 



-iVG(z) 




I gN(a^+g_-V-l) 

1 



1 

-1 



1 

Figure 6. Transformation of parametrix jumps to original jumps 
Also these jumps are described by the same matrices as in the original RHP, indeed 

e'i- 



e 3 



1 1 



1 o\ _ 

,NG{z) J . 2 G ^/ n 74 



;L ^N{9++9^-V-l) 

1 



z G Z^n (/3,oo) 



1 OA -f 



1 
-1 



z eUr\ (a,/3), 



hence the ratio of S{z) with U{z) is locally analytic, i.e. 

S{z) = Nr{z)U{z), \z-p\<r<^. 



(7.7) 



Let us explain the role of the left multiplier Br{z) in the definition ()7.5p . Observe that 
i \ /A, . I —1 i\ I —t I \ I z — a \ - III —I' 



This relation together with the asympotics (|7.3p implies that 



M{z). 



U{z) = M{z) 
1 + 



35 



48(3/2 \di -IJ ' 4608C6/2 \-12i 



-1 12i 



1 / 7(5^2 _ -(7^-2 ^ 5^2) \ 35 / -1 12i 

96(3/2 [i(76~^ + 5^2) -(7^2 - 5^2) J + 4608(^/2 V-12i -1 



+0 



M{z) 



(7.8) 
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as — )■ oo and < n < |z — /3| < r2 < f (so \C\ — )• oo). Since the function ({z) is of 
order N'^/^ on the latter annulus and 5{z) is bounded, equation ()7.8p yields the matching 
relation between the model functions U{z) and M{z), 

B 

U{z) = {1 + o(l))M(z), iV ^ oo, < ri < |z - /3| < r2 < I 

which is crucial for the succesful implementation of the nonlinear steepest descent 
method as we shall see later on. This is the reason for chosing the left multiplier Br{z) 
in (j7.5p in the form (|7.6p . 



8. Construction of a parametrix at the edge point z = a 



This time the construction is similar to the one given in the last subsection. First 
introduce for G C 

and obtain for < arg ^ < vr as ^ — >• oo 



-0-3/4 



MC) 



+o(c-9/2) 



1 i 



-1 6i 



+ 



35 



48C3/2 IJ 4608C6/2 \-l2i 1 



1 12i 



'01 



as well as for vr < arg ^ < 27r 



MC) 



-1 6i 



48C3/2 V 6^ 17 4608C6/2 



+ 



35 



1 12i 
-I2i 1 



Next, instead of ()7.2p . define 



io(C). 



00. 



1 

-1 1 
1 1 
1 



1 
1 1 



arg C G (0, f ), 

arg C e (f ,7r), 
arg C G f ), 
argCG (f ,2^). 



U) 



which solves the RHP of Figure [7] 

• i^^(C) is analytic for C G C\{arg C = 0, f , vr, ^} 
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Figure 7. The model RHP near z = a which can be solved explicitly 
using Airy functions 



• We have the following jumps on the contour depicted in Figure [7] 

arg C = 



ir(c) 




ir(c) 









1 
-1 



1 

1 1 

1 1 

1 



arg C 



3' 3 
arg ^ = TT 



A similar argument as given in the construction of A^^ {(^) implies 

-0-3/4 

Ab.h{C) 



+- 



2^ 
35 



1 i 



I 



-1 Qi 
48(3/2 \ Qi 1 



1 ^^^„9/2^ 



e 3 



C — > 00 



2) 



4608C6/2 \-l2i 

valid in a full neighborhood of infinity. 
Again we use the model function A^^ {(,) in the construction of the parametrix to the 
solution of the S-RHP near z = a. Instead of (j7.4p 



3iV \ 



2/3 



251(2) + V{z) + l + 2TTi sgn(Imz; 



2/3 



z — a < r. 



(8.3) 



This change of the independent variable is locally conformal 

2N ^ 



C(^) 



(z — a) (1 + 0(z — a)) , |z — a| < 



_ {—a)^/|3 — a 

and allows us to define the left parametrix X^z) near z = a by the formula: 

V{z) = S^(z)iy^a3i^^(C(2))e§^'^'^'(")"^ \z - a\ < r 
with the matrix multiplier 

0-3/4 

Xiz)^^ ] , Bi{a) 



B, z) 



-I I 
1 1 



z — a 



-I I 
1 1 



—a 



•4) 
.5) 
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Similar to the situation in the last subsection, V{z) has jumps on the contour depicted 
in Figure [8] which are described by exactly the same jump matrices as in the S'-RHP, 
hence the ratio of parametrix V{z) with S{z) is locally analytic 



a 



(8.6) 

and the left multiplier (jS.Sp in ()8.4p provides us with the following asymptotic matchup 



S{z)=Ni{z)V{z), |z-a|<r<Y 



I gJV(9++9--V-;) 

1 




1 



1 

~1 



iVG{z) 



Figure 8. Transformation of parametrix jumps to original jumps 
between V{z) and M{z) 



Viz) 



M{z) 
1 + 

+o(c 



/ + 



-1 6i 



+ 



35 



48C3/2 \ Qi IJ ' 4608(6/2 \-12i 1 



1 



12i 



5(5- 



i(55-2 + 7^2) 



96(3/2 \^i(5r2 + 752) _(55-2 _ 7^2) J + 4go8^6/2 



35 



1 

-12i 



12i 
1 



-9/2N 



M{z) 



•7) 



as iV 



oo and < ri < \z — a\ < r2 < thus 

^(z) = (/ + o(l))M(z), iV^oo, 0<ri<\z-a\<r2< 



a 



At this point we can use the model functions M{z), U{z) and V{z) to employ the final 
transformation. 

9. Third and final transformation of the RHP 
In this final transformation we put 

(y(z))^\ |z - a| < r, 



R{z) = S{z) { 



(9.1) 



{U{z)) \ \z-l3\<T, 
{M{z)) , \z - a\ > r, \z - (3\ > r 

where < r < min {^j f } is fixed. With C^,/? denoting the clockwise oriented circles 
shown in Figure [9l the ratio-function R{z) solves the following RHP 

• R{z) is analytic for z G C\{Ca,(3 U F U (— oo,a — r) U (/? + r, oo)} with F = 
7+ U U 73" U 74- 
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73 74 

Figure 9. The jump graph for the ratio function R{z) 
For the jumps, along the infinite branches (— oo, a — r) U (/3 + r, oo) 



R+{z) = R^{z)M{z) 







1 



[M{z))-\ 



on the vertical line segment (— ie, ie) 

R+{z) = R^{z)M{z)(^.^^^^ j)(M(z))-\ k 



1,2, 



on the upper lense boundary 7^*" U 7^ respectively lower lense boundary 73 U 74 
R+{z) = R.{z)M{z)I^^JGiz) °)(M(z))-\ ZG71+U72+ 

R+{z) = i?_(2)M(z) (M(z))-\ ze73-U74- 

and on the clockwise oriented circles Ca,i3 

V{z){M{z))~\ \z-a\=r, 
U{z){M{z))~^, \z-f3\=r. 

• In a neigborhood of infinity, we have R{z) — )■ /. 

Here, by construction, the function R{z) has no jumps inside of Cq, and C/3 and across 
the line segment in between. It is important to recall the previously stated behavior of 
the jump matrices as N ^ 00. In fact, on the lense boundaries, the vertical line segment 
{—ie, ie) and the half rays (—00, a — r) U {(3 + r, 00) the stated jump matrices approach 
the identity matrix. Also in virtue of ()7.8p and ()8.7p the same holds true on the circles 
Ca,i3, together, with Gr denoting the jump matrix in the latter ratio-RHP and S/j the 
underlying contour. 



R+{z) = R^{z) 



\\Gr - -^llL2nL«=(Efl) < 



IniV' 



N 



00 



(9.2) 



with a constant c > whose value is not important. The latter estimation enables us to 
solve the ratio-RHP iteratively. 

10. Solution of the RHP for R{z) via iteration 

The stated RHP for the function R{z) 
• R{z) is analytic for z G C\Ilij. 
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• Along the contour depicted in Figure [9] 

R+iz) = R^iz)GRiz), zGJ^R. 

• As z — oo, we have R(z) = I + 0(z~^) . 
is equivalent to the singular integral equation 

R-(z) = I+^ f R^{w){Gr{w)-I)^^ (10.1) 
Ztti J ' w — 

and by standard arguments (see [11]) we know that for sufficiently large N the relevant 
integral operator is contracting and equation (jlO.ip can be solved iteratively in L^(S/j). 
Moreover, its unique solution satisfies 

l|i?--/||L^(E«)<j^, iV^oo. (10.2) 

Observe that for z € CxHr 

R(^z) = I + -^ [ R-{w)(GR{w)-l)dw + 0(z^^), z^oo (10.3) 
zvrz J 



s 



R 



and also as — oo following from (|1U.2|) . (|7.8|) . (|8.7|) and (|9.2|) as well as the previous 
discussion about exponentially small contributions 

ie 

j R-{w){GR{w)-l)dw = j {GR{w)-l)dw+ j {Gr{w)- l)dw + 0{N^^). (10.4) 

Sfl Cq,;9 —ie 

We are now ready to prove the statement on the large A^ asymptotics of /iat given in 
Theorem 1. 

11. Asymptotics of - proof of theorem 1 
We recall the following identities, see p.ip . 

with 

yW(z) = f 7+^ + 0(^-2) j2«-3^ z^oo. 

and trace back the transformations 

Y{z) = y(^)(z) ^ T{z) ^ S{z) ^ R{z). 

Thus 

y/'^)= hm fz(yW(z)z-^'^3-/)) = hm (z{e^''m{z)M{z)e''^3{z)-^,)a,^-Na, 



and since 



giv^g^.;-tja3^-iva3 ^e--r'^3(j_ - -v— "V 3 + 0(^"')), Z ^ OO 
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and similarly 

M(.) = (7 ;) (_\ ::)=/- ^..a + 0(.-), . ^ oo 

we have in combination with (jl0.3p and (jl0.4p as — )• oo 

4 4 zvr J 



is 



= -Ar(/3 + a)^-(/3-a)^ + ^ y {Gr{w) - l)dw + ^ J {Gr{w) - l)dw 

+0(iV-2) (11.1) 

We start computing the contributions arising from the circles Ca^p. From (j7.8p as — > 

oo 

y - I)d^ = ^ y (^,(75-2 + 5^2) _\75-2 _ 5^2^) J ^57^ + 0{N ), 

and since 

. y/3 , 



for z : |z — /3| < r, we obtain via residue theorem 

y = " 7 = " 4iV V' " 5/3(/3 - a) 

Combined together as A'" — )• oo 

(G/jH (11.2) 

27ri / 3/52-2a/3-a2 ^(11/3^ - 12a/5 + a2)\ . 2^ 
192Af(/5-a) V^(ll/32-12a/3 + a2) -(3/52 _ 2a/3 - a^) ^ + ^^^^ ^ 



For the integral over Cq we follow the same strategy. First from ()8.7p as — )• 00 
/■^^ . . ^ 5(^-2-7^2 «(5r2 + 752)\ 

and with 

iTr/^ 2Ar \2/3^ 5a/3 - 3a2 - 2/32 s2n^^^^ 

= ' ( J - + 10(-a)(/3-a)2 - ^) + 0((. - a)2) 

for z : |z — a| < r, we deduce from residue theorem 

5-2 (-u;) , (-a)i/ 5a/3-3a2 -2/32\ /■ ^2^^) 

^^du; = 2..^ ^1 - J , y ^^^^ = 
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SO together 

j {Gb{w) - I)dw (11.3) 



192A^(/3 - a) - 12a/3 + /S^) Sa^ - 2a/3 - /S^ 

Adding (fTL2]l and (fTL3]) we have thus 

Cq,;3 

Let us now move on to the evaluation of the integral 



l£ 



{GRiw)-l)dw = _,(52 _ J 









+ J 4 [{5 + 6-')'' -i(52 - r2)J ^^^-^^ 

—ie 

As we see this evaluation requires the asymptotical computation of the integrals 

ie 

ji{w)f{w)dw / j2{w)f{w)dw (11.6) 



for a function /(tf) which is analytic on {—ie, 0) U (0, ie). Consider the first integral, we 
have from (|5^ 



ji{w)f{w)dw = -2e-^'^2{0) / g-7V(-|2,in3,+s/hi(y)+ft2(y)-h2(0)) sm{Ny)f{iy)dy 





where 



, / N 4 ^/l3{iy -a) + -a{^ - iy) yjiy-a + i^J (5 - iy 

hi{y) = -\n 7^= ' /i2(y) = 41n 7^= • 

vr p — a VP — a 

Let us perform the change of variables 

u = u{y) = --y\ny + yhi{y) + h2{y)-h2{Q), yG[0,e). (11.7) 
vr 

We have ii(0) = and would now like to express y as a function of u. To this end 
introduce v = u/y, thus from ()11.7p 

''.\nu+-\nV + hi{-)+-(h2{-)-h2{<d)\ (11.8) 



vr vr W/ u 

and we are going to solve this equation for v by iteration 

fo = -— Inii, Vn+i = -— Inu + — lnti„ + /iif— ) +—{ h2( — ] - /i2(0) 
vr vr vr \v„ u \ \v„ 
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First 



Vl 



n 



\nu+^ln(^ - ;^lnn) + hi{0) + h'^iO) + O 



u 



and secondly 

f 2 = - 



2 2/2 

-lnti + -ln Inn) + hi{Q) + h^i^) + O 

TT vr \ vr 



In It 



ln(— Inn) 
— Inn 



This asymptotic behavior persists for subsequent v^s, we have a solution to equation 
(fTLH]) of the form 



2 2/2 

v = lnn + -ln Inn) + hi{0) + h^iO) + O 



IT 



vr 



vr 



ln(— Inn) 
— Inn 



as n — 7- 0. This in turn implies for the solution y = y{u) that as n — t- 0, 



n 



y 



--Inn 

as well as for its derivative 



^ In(-flnn) + Ki^) ^ ^// ln(-lnn) Y 

— In u — - In u V I — In zi y 



du 



2 2 

- In y - - + /ii (y) + yh[ {y) + /I'a (y) 

TT TT 



(11.9) 



(11.10) 



1 



In n 



^ In(-flnn) ^ I - ^{hm + Ki^)) , p// ln(-lnn 



Inn 



Inn 



Inn 



At this time we go back to the given integral 



j h{w)f{w)dw 





u{e) 

_2^-Nh,{0) f e~^^sm{Nyiu))f{iy{u))^du 



and introduce 



A{N) = 4 / e-^''-^y^''y+y^^^y^+''^^y^-^^^^^hm{Ny)dy. 



(11.11) 



Proposition 1. As N ^ oo 

1 



A{N) 



Ar(lniV)2 



^ _ 2 In In TV 3 - 27 - 2 In (f) - 7r(/ii(0) + h'^jO)) 
IniV hT/V 



+0 



InlniV 
InA^ 
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Proof. We use the change of variables (jll.9p . (jll.lOp as weh as the substitution s = Nu, 



Nu{e) 

1 r se-' 

N J (lniV-lns)2 




21n(^lniV- ^Ins) 
In — In 5 



1-^(^1(0) + ^2(0)) ln(flniV-fln.) / , w ..^ 
lniV_ln. (lniV-ln.)2 V ' + 

In(lniV-lns) 



-O 



In — In s 



ds 



Since u{£) can be chosen arbitrarily small, we can expand the integrand in reciprocal 
powers of In A^. This gives as — )• oo 



A{N) 



1 



Af(lniV)2 



, , 21n(^lniV) , , 



InN 



I l-7r{h{0)+h',{0)) ^^ ^ 2h{N) ^ ^ / /In In iV 



IniV 



IniV 



IniV 



with 



Nu{£) Nu{e) 

Ii{N)= j e-'sds, l2{N)= j e-'shisds. 



Up to an exponentially small error, we have Ii{N) = 1 and l2{N) = 1 — 7, where 7 
denotes Euler's constant, thus 



AiN) 



1 



A^(lniV)2 



_ 2 In In TV 3 - 27 - 2 In (f) - 7r(/ii(0) + h'^jO)) 
IniV ^ hTiV 



+0 



In IniV 
IniV 



, iV 00. 



(11.12) 



□ 



Back to the integral under consideration, we notice from pi.9p 



/{■'K^))-wo)=oy 

hence via Proposition [1] 



InN 



, iV — ;> 00, 



ie 

I n{w)f{w)dw = -ye-^'^^°V+(o)A(iV) + o(^^^^l^) (11.13) 
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with 



hi{0) = - In , /i2(0)=41n = 4i arg , 

vr VP — a -y/p — a vP — a 

The second integral in pi.6|) can be treated in a similar way. Indeed we have 

£ 

J 32{w)f{w)dw = -2e-^'^^(°) j e-^(-|^'''^^'+^''^3(^;)+/^4(^/)-/^4(0)) ^:,^^Ny)f{-iy)dy 

-IE 

with 



, ,x 4 ^/f3{-iy -a) + y^-a{P + iy) V"^?/ - « + W/^ + 

/i3(y) = -ln 7^= ' /i4(y) = -41n 7^= 

vr VP — a VP ~ 

and we deduce as iV — oo 



— ie 

with 



/i4(0) = — 4i arg 



\//3 - a 

Adding (|11.13p and ()11.14p . we end up with 

ie 

1 j2M/Md^^ = -y(e-^''^W/+(0) + e^'^^W/„(0))^(iV) 

-ie 

The latter expansion enables us now to evaluate (jll.Sp . Since 

1/4 / \ 1/4 



6 4 



we obtain 

{GR{w)-l)dw = -^A{N) X 



+0 



_2^((^)i/2 + (^)"i/2) sin^^ -4cos<^;v - 2((^)i/2 - 
4cos(^jv - 2((^)i/2 - gi^^^ 2f((^)V2 + (^)-V2) sin^^ 

1 



iV2(lniV)3 



EXACT SOLUTION OF THE SIX- VERTEX MODEL 



31 



as — 7- oo with 



All together from (fTTT]) 



ipN = arg- 



-a 



a 



^^,3 (7V)^,3__ 1 / /3 + a 3i(/3-a) 



16 



A{N) X 



1^4cosv9;v - 2((^)i/2 - sin<^^ 2i((^)V2 + (^)-V2) ^in^^ 

and in particular 



(N) 



12 



„Nl 



+2 



a\i/2 /-a\-i/2 



/3 



/3 



4 cos AT 

}sin(/p^^ +0{N~^) 



The latter expansion allows us to deduce the asymptotics of the normalization constants 
Hn as —7- 00. Since 



2N+1 fvW 



12 



one obtains 
hN = 



1 TTAiN) 
1 + TTT + 



X ( 4cOS (/3Ar + 2 



{(f) 



AN 4(/3-a) 
a\i/2 /-a\-i/2 



I sin LpN^ 



+ 0{N 



-2> 



iV ^ 00 



or in other words, recalling the definitions of a, (3 and I and the identity 

1 



(PN = ANavg ^== — = 4A^arctan 



a 



we deduce 



2Vcosf^y 16^ 



tanf (1 - x) 
1 , (-1)^ 



7riV(l 



^ + 4iV+ 2 



cos (7rx( N + -] ]A{N) 



+0(iV"2) 



(11.16) 
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Here the stated expansion as — )• oo is uniform on any compact subset of the set (jl.lSp . 
Furthermore by Stirhng's approximation 

(11.17) 



one obtains 



= r e 6^ 

(A^!)2 V2cos^y 



with 



= ^ + ^^cos(^x(iV + i))A(iV) + o(^^-2) 

1 (-l)A^cos(7rx(iV+i)) r 
AN 2iV(lniV)2 \ 

thus proving Theorem [TJ 



21nhiAr l-27-41n2-21n(cosf^) 
InN MiV 



12. ASYMPTOTICS OF THE PARTITION FUNCTION - PROOF OF THEOREM [2] 

We go back to (fLT2|) 

and derive in the given situation from ()11.17p 

Zj, = CF^\^^^->^-^\ F = 4i^ (12.2) 

^ COS 2 

with a positive, independent, constant C. Applying now Euler's summation formula, 
we have as — )• oo 

AT-l AT-l 

^ (^'^^ " 6fc ) = 12 ^ + "° + 2 ^ ("""^ + 2 ) ) + ^ (^"'^ ' 

fc=i fc=i 

with an independent term cq. The sum can be further estimated using summation by 
parts: 

7V-1 

'^{-1)'' COS {7rx( k + -))A(k) = A{N-l)SiN-l) (12.3) 



k=l 



N-2 



J2 {A{k + l)-A{k))S{k) 



k=l 



with 



S{k) = X:(-l)' cos {.x{l + l)) = (-l)^^o^^^(fe + l))-cos^^ _ 

1=1 2 ^ 
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Following the notation of Proposition [H as — t- oo 

u(e) 



hence 



u{e) 

A{N + l)-AiN) = ^/e-^"(e-"-l)-n(A^|/(n))|d. + 0(^ 



NY 



^ >^(logiV)^ J' ""^"^ 



and therefore the series 

oo 

Y^{A{k + l)-A{k))S{k) 

k=l 

is absolutely and uniformly convergent on any compact subset of the set (jl.lSp . Back 
to (|1'2.3|) using Proposition [T] 

N-l 

Y,{-^?cos {7:x{k + -))A{k) = A{N -l)S{N -l) + Co 



oo 



+ {Mk + 1) - A{k))S{k) = Co + o(^^^j^) , iV ^ oo 

with an N independent term Cq, thus 

Af-l 

= U^''^ + ^ + 0{N-'), N^oo. (12.4) 

k=l 

Back to I^TT2\\ . we get as iV oo 

72 „ 1 



Zn = CF^ N-2(l + 0(N-^)]. (12.5) 



As mentioned before, the x dependency of C will be derived from the Toda equation 
(fL9]) . It imphes 

d'^ , TN+lTN-l hN , . 

dx^ hN-i 
and we now use equation (jll.l6p 

" 21nlniV 



7r/27riV\2^+ie-2" 



2 V cos ; 16*^ 



. , 1, (-l)^cos(^x(iV + l)) f 
4A^ 2Af(lniV)2 \ 



co(x)^^//lnln^yxl ci^^^ 



InTV VV IniV J J ] \N'^{lnN)'^ 
uniformly on any compact subset of (jl.lSp . Here 

/ TTX 

co(x) = 1 - 27- 41n2 - 21n COS 



InA^ 
iV ^ 00 (12.7) 



2 
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has already been computed, but not ci{x). Substituting (jl2.7p into ()12.6p . we get 



hN-] 



TT 



2cos^ 



'N{N -1) 



1 \N 



InN 



InN 



In In N 



1 + 

2 



{-ly^ COSTTXN COS^ 

Af(lniV)2 



4iV2 



+ 



Af2(lnA^)2 



IniV 
N ^ oo. 



Since 



1 + 



X \ 
N 



) 



N ^/ x^ 3x^ + 8x3 ^.,,_3 



we can simplify ()12.8p further 



hN-i 



irN 



2cos^ 



.0{ 



In In TV 



IniV 

At this point we notice that 



( — 1)" cos vrxA^ cos ^ 
iV(lnA^)2 

' +0 



12Af2 



vr 



2cosf^ 



In cos 



iV2(inAr)5 



TTX 



(12.8) 

, N oo 

21nlniV co(x) 
IniV ^ InN 

N -^oo. (12.9) 



as well as 



vr 



2cos^ 



TTX 

COS TTXjyj COS 

2 



which implies with (jl2.9p 

hN _ 



■KX 



cos vrxA^ 
4A^2 cos 



In cos — I (N^ - — ] + 



(InNy 



and the error term in this identity is uniform for any x chosen from a compact subset 
of (fLT5]) . Thus from (fTZTjl 

and by integration with respect to x 



1 



1 



In H In cos 

12 /12 V 2 



TTX 



+ 



1 



(lnA^)2 



InZ^ = In ""^^ fj + 1 IniV + Incos ^ + di(A)x + do(A) 

1^ COS Q -L^ ^ 



+0 



1 



(12.10) 



,(lnA)2^ 

where do and di in general depend on N, but not on x. Substituting ()12.5p into ()12.10p . 
we obtain that 

1 



1 TTX 

InC = — Incos — + di{N)x + d2{N) + O 



(12.11) 
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and C does not depend on N . For any xi,X2 G (—1) 1) this implies 

In C(xi) — In C(x2) = — In cos Incos — ^ + (ii(A^)(xi — X2) 

12 2 12 2 



+0 



1 



(IniV)V 



i.e. the limit 



exists, hence also the limit 



lim di{N) = Ki 



lim d2{N) = K2 
N—>oo 



exists. By taking the limit — )• 00 in ()12.1ip . we obtain that 

1 TTX 

In C = — In cos — + kix + K2 

thus proving 

Proposition 2. The constant C in the asymptotic formula (|12.5p has the form 

C= ( cos^V'e"^"+"^ 



We will now show that, in fact, ki = 0. To this end recall the initial Izergin-Korepin 
formula (|1.7p and (jl.Sp . in our situation 

(1 _ 2;2^Af2 / \N ^ 

= TTTlvrTTTTTT^JV' = det , ■ , ip{x) 



It shows, that as a function of x is even, hence in the notation of Proposition [2l we 
conclude 

Hi = 0, 

thus proving Theorem [2l 

13. Phase transition - proof of Theorem [3] 

We start with change of variables (|1.3ip . It implies in the disordered phase region, 
via (fT3|l . 

sin(7 — t) 1 — X sin(7 + t) 1 + x 

sin(27) " ^ ^' sin(27) " ^ ^' 

hence for y > 0, 

/ y(l + y) . 
sm7 = 2W ^-^^^^^2_^2 ' smt = 2;sm7. (13.1) 

On the other hand, in the antiferroelectric phase region, via (|1.4p . 

sinh(7 — t) 1 — X sinh(7 + 1) 1 + x 

sinh(27) ~ 2 ^' sinh(27) ~ 2 ^' 
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hence for y < 0, 



smh7 = 24/- ^ n 1 sinht = X sinh7. (13.2) 

y (1 -|- 2yY — 

The functions sm z and sinh z are both entire, satisfying the usual relations, 

sin(— z) = — sin z, sinh(— z) = — sinhz, sinhz = — isin(iz). 

This implies that the inverse functions, arcsinz and arcsinh z, are analytic at z = 0, and 
they satisfy the relations, 

arcsin(— z) = — arcsinz, arcsinh(— z) = —arcsinh z, arcsinhz = — z arcsin(zz). 

Let us make the change of variable y = /c^ in the disordered phase and y = — /c^ in the 
antiferroelectric phase, where A; > 0. Then p3.ip implies that 



7 = /g(x, k) = arcsin ( 2k A ^ ^^^^^ _ J \, (13.3) 



while from (jl3.2p . 



7 = /af(2;, k) = arcsinh \^2k^ j . (13.4) 

Here the both functions, foix, k) and fApix, k), are analytic at A; = for any x G (—1,1), 
satisfying the relations, 

fD{x,-k) = -fD{x,k), fAFix,-k) = -fAF{x,k), fAF{x,k) = -ifD{x,ik). 

Now, from (|13.1|) and (|13.2|) we have that 

t = gD{x,k) = arcsin (^2xsm^^ = arcsin (22; sin /£)(x. A;)) (13.5) 

i = 9AF{x,k) = arcsinh (2rc sinh 7) = arcsinh(2x sinh/^i?(2;, /c)) . (13.6) 

Observe that for any x E (—1,1), the functions go{x^k),gAF{x^k) are both analytic at 
A; = 0, satisfying the relations, 

gD{x,-k) = -g{x,k), gAF{x,-k) = -gAF{x,k), gAF{x, k) = -igD{x,ik). 

The four equations, (|13.3p . (|13.4p . ()13.5p . and (|13.6p . prove stated behavior ()1.34p of 7 
and i as y — )■ 0. 

Next we go back to ()1.23p and ()1.30p and employ change of variables (|1.3ip . 
1^ / 1 2\ ™^ sin(27) 

FD{x,k) = ^- ^, 7 = /£,(x,A:), t = gD{x,k) 

27 cos 

rfeg/ ,2\ vrao smh(27j 

F^p,{x,k ) = ^, j = fAF{x,k), t = gAF{x,k). 

27 cos 

We first notice that the functions 

sin(27) sin2/z)(x. A;) sinh27 sinh2/AF(a^, A;) 

hnix, k) = — = — — , hAF{x, k) = — = — - — - — -— 

27 2fD[x,k) 27 2fAF{x,k) 
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are both analytic at A; = and 

hD{—k,x) = hD{x,k), hAF{—k,x) = hAFix,k), hAF{x,k) = h£){ik,x). (13.8) 
Similarly, the ratios 

t _ gpjx, k) t _ gAF{x,k) 

7 fD{x,k)' 7 fAFix,k) 
in the corresponding phase regions, are analytic at A; = satisfying the same relations 
as in (HMD, 

gpjx.-k) ^ gpjx.k) gAFjx, -k) ^ gAF{x,k) gAF{x,k) ^ gpjxjk) 
fD{x,-k) fD{x,k)' fAF{x,-k) fAF{x,k)' fAF{x,k) fD{x,ik)' 
But this shows that the functions 

_ sin(27) _ sinh(27) 
"^oix, k) = ^, rAF{x, k) = ^ 

27 cos 1^ 27 COS ^ 

are analytic at /c = and they satisfy again the same relations, 

rD{x,-k) = rD{x,k), rAF{x,-k) = rAF{x,k), rAF{x,k)=rD{x,ik), 

so in a neighborhood of A; = we have the following convergent Taylor expansions: 

00 00 
rD{x,k) = ^rj(x)A;2j, rAF{x,k) = J^i-'^Y rj{x)k^^ . (13.10) 
j=o j=o 
Substituting back y, and noticing that 

irab /I \2 
^=(2 + 
is clearly analytic at ?/ = 0, we obtain 



Fnix,y) = r^'^{x,y) = Y,fji^)y', x G (-1,1), (13.11) 

with 



, 7r(l-2;2) ^ 7r(7rx'^sin^ -vrxsin^ + 8cos^ 
^'^^^ = T^^ ^^(") = 12cos2^ ^ 



2 J-^^wo 2 

and the stated series in (jl3.1ip is convergent for small y. 
Finally, to prove (|1.37|) . observe that by (jl.28p . 



F^^^(7,t) = 0(g2) = o(e-4) 
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